Instanton approach to the Langevin motion of a particle in a random potential 
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We develop an instanton approach to the non-equilibrium dynamics in one-dimensional random 
environments. The long time behavior is controlled by rare fluctuations of the disorder potential 
and, accordingly, by the tail of the distribution function for the time a particle needs to propagate 
along the system (the delay time). The proposed method allows us to find the tail of the delay time 
distribution function and delay time moments, providing thus an exact description of the long-time 
dynamics. We analyze arbitrary environments covering difi^erent types of glassy dynamics: dynamics 
in a short-range random field, creep, and Sinai's motion. 



Introduction-Owe dimensional driven dynamics in a 
random environment attracts a great deal of current at- 
tention. The motivation of the interest is two-fold: First, 
the propagation of a particle through a ID random field 
has become a paradigm for a general out-of-equilibrium 
stochastic processes in random systems capturing all ef- 
fects of pglassy dynamics including aging and memory 
effectsETQ. Second, a particle moving in a ID random po- 
tential models straightforwardly a variety of physical sys- 
tems ranging from dislocations and charge density waves 
in solids, spin-chain dynamics and doraain. growth, to 
protein molecules and bacterial coloniead^ErQ^El. The at- 
traction of the ID models is thus that while allowing for 
an analytical treatment (and often even for a full analyt- 
ical solution) they also offer an insight into generic ba- 
sic properties of the wealth of glassy out-of -equilibrium 
systems. Indeed, even the simplest ID models with the 
Gaussian correlated potential v{x), with {v{x)) = 0, 
{{v{x) — 1;(0))^) = K\xy exhihit a striking generality and 
diversity of glassy behaviorsB; 1 

The approach employed iroB enables the exact deriva- 
tion of the particle velocity, but it does not allow for 
a complete dynamic description of the system - for ex- 
ample, in terms of velocity cumulants and/or correlation 
functions. Taking a kinetic view on the problem one can 
characterize the dynamic properties of a random system 
by the probability distribution for the particle velocity 
or, equivalently, by, the probability distribution for the 
particle delay time r. The latter is defined as an average 
time that a driven particle spends to propagate through 
the sample, and the corresponding distribution function 
P(r) characterizes completely the transport properties of 
the ID system involved. 

pXhe exact expressions for a particle velocity derived 
inl9fl indicate that the system dynamics is governed by 
rare fluctuations of the random potential, i.e. by the tail 
of the delay time distribution function. Thus the instan- 
ton solution for the Langevin equation in the presence of 
an external force offers a most adequate description of the 
long time glassy dynamics in an arbitrary ID random en- 
vironment. In this Letter we apply the instanton method 
developed earlier inQ to find the asymptotic behavior of 



the delay time distribution function P{t) at large t. In 
general, the idea of the instanton method is to pick up 
the largest contribution to the functional integral coming 
from an optimal configuration instead of implementing 
the complete integration. In our case of low enough tem- 
peratures, the delay time is mainly determined by the 
transition over the largest barrier in the system. Since 
this time is exponentially large, one expects that the de- 
lay time averaged over the random potential (or any of 
its higher moments) is determined by the saddle point 
trajectory of some effective action controlling this expo- 
nential behavior, and, thus, it can be most appropriately 
found by the instanton method indeed. Further, knowing 
all the moments of the delay time one can reconstruct the 
large r asymptotics of the distribution function P(r). 

The distribution function P[t) depends strongly on 
the specific form of the correlation function of the random 
potential u{x). We restrict ourselves to the Gaussian 
disorder with the correlation function 

{v[xi)v{x2))d = u{xi - X2), {v{x))^Q. (1) 

We use also the function 

K{x) = {{v{Q) - v{x)f)d = 2{u{Q) - u{x)), (2) 

which is more convenient in the case of the long-range 
correlated potential. In the case of the short-range cor- 
related potential, when u{x) monotonically decreases to 
zero as a function of |a;|, we find that the delay time dis- 
tribution is log-normal 



InP(y) 



'4w(0) 



(3) 



where y = Inr. The distribution function -P(F) is de- 
fined as the distribution function of In r, being therefore 
related to the distribution function P(t) by 



P(r) 
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(4) 



In the case of the long-range correlated potential K{x) = 
K x'' , with K > and < 7 < 1, the distribution function 
essentially depends on the applied field E: 
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And, finally, for the so-called Sinai modelS K{x) = kx 
we obtain: 



P(r) 



. -2TE/K 



(6) 



which agrees with the result found earlier in Ref.B (see 
alsoEl) 

The model and results- Now let us turn to the detailed 
calculations. The dynamics of a classical particle in a 
random potential is described by the Langevin equation 



T-^dtxit) = -l3d,v{x) + l3E + ^{t), 



(7) 



where (3 is the inverse temperature, F is the inverse relax- 
ation time, v{x) is the random potential, E is the applied 
uniform field and ^(t) is the Langevin thermal noise that 
models the thermal environment. 



(e(ii)C(i2))T = 2r-i(5(ii-t2). 



(8) 



To distinguish between the average over the thermal 
noise and disorder averaging we denote the former by 
the subscript T. The inverse relaxation time F will be 
set henceforth to 1 for convenience. Making use of the 
standard approach (see for example Ref.tj) one can write 
the Lagrangian corresponding to Eq. (|^) in a form: 

jC = -x^ - ix {dtx -I3E + f3d^v{x)) + ^v"{x). (9) 

Further we will set (3 = 1 measuring energies in the units 
of temperature. The probability for a particle to go from 
the point xi to X2 is given by the functional integral 



P{x,,X2)= I D[x{t)]D[x{r)\<i^''"'. 



(10) 



and the time a particle spends moving from xi to X2 (the 
delay time) is, correspondingly. 



t{xi,X2) ^ P ^{xi,X2). 



(11) 



If the temperature is lower than a typical barrier in the 
system, one can use the saddle point approximation to 
find P 



P 



A 



dt Cs 



(12) 



where Cg.p. is the value of the Lagrangian (g) at the sad- 
dle point trajectory. The delay time t{xi,X2) averaged 
over the disorder is then 



with the effective action 



(13) 



A[v, = - 2 y dxidx2 v{xi)f{xi,X2)v{x2) - J dt C, 

(14) 

where / is the inverse correlation function (|^): 
/ dx f{xi, x) u(x, X2) = (5(xi — X2). Taking the variational 
derivatives of the action (llj) with respect to i, x, v we 
get the saddle point equations 



-2x + dtx- E + d^ vix) = 0, 
dtX — xv"{x) = 0, 

v{x) = — / dtx{t)u'{x — x{t)), 



(15) 
(16) 

(17) 



where the variable x was redefined ix —>■ x. When deriv- 
ing Eqs.( p^ , p^JT7| ) the last term in Eq(||) was neglected 
because it gives the 0(1) contribution to the action while 
the whole iustanton action is much larger than one. Fol- 
lowing Ref.Ll, in order to find an instanton solution we 
have to choose x = dtX. Indeed, this allows to reduce 
Eqs.(E5|,E6|) to a single equation: 



dtx = —E + dxv{x). 



(18) 



Assuming that the instanton solution x = dtX exists in 
the interval {rji,ri2) from Eq.(p^ we find 



v{x) = u{x — 772) — u{x — 771). 



(19) 



The interval (ryi, 772), where the instanton solution exists, 
obviously must lie within the interval (xi, X2)', we will 
see later, however, that these intervals do not necessarily 
coincide. Outside the interval (771, 772) one should take 
the normal solution in the form 

£ = 0, dtx = E~dxv{x). (20) 

The action (|lj) corresponding to the solution ( |l^ ) is 

1, 



(21) 



Since within the interval (771,772) x = dtx and outside 
this interval x = 0, assuming that x is continuous, we 
conclude that on the boundaries of the interval (771,772) 
X = dtX = 0. Therefore from Eqs.( p^ , p^ ) we get the 
equation 



E 



-u'{r]2 - 771) 



(22) 



which defines the distance 772 — 771. Using Eq.(p^ the 
action ( ^l|) can be written as 



A = it(0) - u(77) - Er], 



(23) 



where 77 = 772 — t^i . Note that Eq. (g2|) corresponds to the 
extremum of the action ( |2^ ) with respect to 77, so the 
final answer for the action can be conveniently presented 
as the maximal value of the action 
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(24) 



where K{x) is given by 

The formula (^) allows to find the average delay time 
for an arbitrary disorder correlation function u{x) and 
it's concrete realization follows from the specific form of 
the correlation function u{x). 

(i) Short range potential. Let us begin with the sim- 
plest case when u{x) monotonically decreases to zero as 
a function of The corresponding form of the effective 
potential v{x) is shown on Fig. 1. From Eq.(p^ it fol- 
lows that 7] oo when E ^ 0, thus if E is low enough, 
the action becomes 



A = u{0). 



(25) 



Tl2 



^1 



FIG. 1. The instanton solution for the potential v{x) in 
the case of short range potential. 

(ii) Potential with long-range correlations. Let us con- 
sider the potential of the form K{x) — kx"^ , with k > 
and < 7 < 1. In this case the expression for the action 
( p^ ) still holds; however, minimizing (^4|) with respect to 
77 we arrive at the action diverging at low fields: 



A 



(26) 



in a|-full agreement with the exact solution obtained in 
Ref.a 

(iii) 'Extremely correlated' disorder, 7 = 1: This is 
the well-known Sinai model K{x) = kx. Now Eq.(|2^) 
is either never satisfied or it is satisfied identically for a 
special value of the applied field E = E^ = k/2. If E < 
Eq, the instanton solution exists, but the initial and final 
points of the instanton motion 771 and 772 coincide with 
points xi and X2 respectively, and the action becomes 



A = {k/2-E) {x2-xi). 



(27) 



It means that the average time always essentially (expo- 
nentially) depends on the boundary conditions. In the 
case E > Eq the instanton solution does not exist. 



Distribution function of delay time. The method de- 
scribed above can be simply generalized for a calculation 
of higher moments of t{xi, X2) 

Tn{xi,X2) ^ {T''{xi,X2))d- (28) 

The action A„ that determines the n*'' moment 

(T"(a:i,a;2))d = e^"['''-] (29) 

is given by 

A„[v,x] = J dxidx2 v{xi)f{xi,X2)v{x2) - n JdtC. 

(30) 

From this equation one can see that An[u{x),x] = 
n A[nu(x), x], therefore using Eq.(|4|) we get the expres- 
sion for the action An 

An=n'^ K{r])/2~nE'q, (31) 

that must be maximized with respect to rj. This gives 

nK'{ri)/2 - E = 0. (32) 

Knowing all moments of r one can find the distribution 
function P{t). Indeed, in terms of the distribution hmc- 
tion P(r), the moments t„ are defined as 



= / P(t) t" dT = / P{Y) e"^ dY. 



(33) 



Since we use the instanton method to find r„, with the 
same accuracy we can use the saddle point approximation 
in calculation of the integral over Y in Eq. ( |33|) . Compar- 
ing Eqs.( |3l|j3^ ) we get 

\hiP{Y)+nY^^ = \v?K{r])/2- Enr^^ , (34) 

where { }y and { }^ mean taking extrema with respect 
to Y and 77 respectively. Differentiating ( ^4| ) with respect 
to n we get 



Y + Er] 

Using Eq.(Q) we find the distribution function 

inP(r) = -li^l±^, 

^ ' 2 K{ri) ' 
where 77 is defined by the equation 

i^'(77) {Y + Erf}-2E K{r]) = 0, 



(35) 



(36) 



(37) 



following from Eqs.( |3^ , p5| ). Interestingly, this equation 
also follows form the extremum of the logarithm of the 
distribution function (pq) with respect to 77. 
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The distribution function (jS^) is our main general re- 
sult. Now we shall analyze it's behavior for different cases 
introduced above: 

(i) In the case of the short-range correlated potential, 
assuming that the applied field is not too strong, from 
Eq.(^6|) we obtain the log- normal distribution 



InP(r) = - 



y2 



4u(0)' 



(38) 



(ii) In case of the potential K{x) — kx"', 0<7<1 
from Eq.(|37|) we get 



and the distribution function (36f) becomes 
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(39) 



(40) 



(iii) Although in the Sinai case (7 = 1) the moments 
of T are not defined in the limit of a large system, the 
Eq.(^0[) has no singularities when 7 = 1. Therefore the 
distribution function in the Sinai case is given by 



Pit) = t 



(41) 



which agrees with the earlier result] 

In conclusion, we have found the asymptotic behav- 
ior of the delay time distribution function for a general 
problem of the Langevin motion of a particle in a one- 
dimensional random potential. The application of our 
procedure to the Sinai model recovers the earlier results 
verifying our approach. The developed instanton method 
allows to derive the distribution function for an arbitrary 
ID random potential, and can serve as an initial step to- 
wards a quantitative study of the glassy dynamics of the 
general multidimensional systems. The method proposed 
can also be generalized for calculation of other correla- 
tion functions which are determined by the contribution 
from the largest barrier in the system. 
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